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We study the current controlled modulation of a nano-contact spin torque oscillator. Three 
principally different cases of frequency non-linearity (d 2 f jdl\ c being zero, positive, and negative) 
are investigated. Standard non-linear frequency modulation theory is able to accurately describe 
the frequency shifts during modulation. However, the power of the modulated sidebands only 
agrees with calculations based on a recent theory of combined non-linear frequency and amplitude 
modulation. 
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Spin-torque oscillators (STO) offer a combination of 
attractive properties such as ultra wide band frequency 
operation,^* 2 - extremely small footprint (without any need 
for large inductors), and easy integration using well- 
established magnctoresistivc random access memory pro- 
cesses. The basic principle of a spin-torque oscillator is 
based on the transfer of angular momentum from a spin- 
polarized current to the local magnetization^^ The effect 
usually occurs in a nanoscale device where a large cur- 
rent density (~ 10 s A/cm 2 ) can drive the precession of 
the free layer magnetization at GHz frequencies,— £ thus 
acting as a nanoscale oscillator. Effective modulation 
of the microwave signal generated from STOs is required 
for communication applications. However, both the STO 
frequency and amplitude are typically non-linear func- 
tions of the drive current. This non-linearity is related 
to a change in the precession angle with the increase 
in the current magnitudei 7 - - — Experiments have shown 
other sources of non-linearities such as temperature^ and 
dynamic-mode hopping.— ~— The wide range of possible 
sources of non-linear behavior is likely to render the fre- 
quency modulation of STOs highly non-trivial. 

Despite the rapidly growing literature on the many dif- 
ferent aspects of STOs, experimental studies of frequency 
modulation are still limited to a single work by Pufall et 
al^ They observed both unequal sideband amplitudes 
and a shift of the carrier frequency with modulation 
amplitude, which they ascribed to non-linear frequency 
modulation (NFM). While linear frequency- modulation 
(LFM) theory assumes that the instantaneous frequency 
of the modulated signal is linearly proportional to the 
modulating signal^ NFM theory takes into account the 
non-linear change in the intrinsic operating frequency 



during modulation. Pufall et al. calculated the ob- 
served sideband amplitudes using NFM theory and found 
a rather large (about 50%) discrepancy between their cal- 
culated and experimentally observed sidebands, which 
they argued might be due to amplitude modulation or 
other non- linear properties of the STO. 

In this work we study the frequency and amplitude 
modulation of a nano-contact STO for various amounts 
of frequency non-linearity. The frequency non-linearity 
is described by the second derivative of the frequency, 
/, with respect to the dc bias current, Id c , d 2 f/dlj c . 
Three different cases of frequency non-linearity (d 2 f/dl% c 
being zero, positive, and negative) are investigated. As 
expected from NFM theory, the carrier and its associ- 
ated sidebands exhibit a change in frequency under mod- 
ulation, which can be directly calculated from the ex- 
perimentally determined non-linear properties of the fre- 
quency of the free-running STO. However, the power of 
the modulated sidebands is only poorly reproduced us- 
ing NFM theory and we show that it is essential to con- 
sider amplitude modulation in order to reach any quan- 
titative agreement. Using a recently proposed theory of 
combined nonlinear frequency and amplitude modulation 
(NFAM)^ we are able to show remarkable agreement 
between our experimental data and calculations, which 
involve no adjustable parameters. Despite the complex 
phenomena involved in the STO non-linearities, we show 
that modulation of these devices is highly predictable. 

The nano-contact metallic-based STOs studied in this 
work have been described in detail in Refii. Us- 
ing e-beam lithography, a circular Al nano-contact 
with nominal diameter of 130 nm is fabricated 
through a Si02 insulating layer, onto a 8x26 fim 2 
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FIG. 1: Current dependence of the free running STO: (a) fre- 
quency, /, and its second derivative, d 2 f /dl^ c ; (b) integrated 
power, both measured in a magnetic field of H = fO kOe, 
applied at 70 ° to the film plane. Dotted lines indicate the 
three different operating points (28, 31, and 38 mA) used to 
compare three principally different cases of frequency non- 
linearity, corresponding to d 2 f/dl% c being zero, positive, and 
negative, respectively. Inset in (b) shows the measured S- 
parameter, Su at the STO. 



pseudo-spin-valve mesa with the following layer struc- 
ture: Si/Si0 2 /Cu(25 nm)/Co 81 Fci 9 (20 nm)/Cu(6 
nm)/Ni 80 Fe 20 (4.5 nm)/Cu(3 nm)/Pd(2 nm). While all 
data presented here has been taken on a single device, 
similar behavior has been observed in several other de- 
vices of the same size. 

The low frequency (100 MHz) modulating current is 
injected from an RF source to the STO via a circu- 
lator. The dc bias current is fed to the device by a 
precision current source (Keithley 6221) through a dc- 
40 GHz bias tee connected in parallel with the transmis- 
sion line. The signal is then amplified using a broadband 
16-40 GHz, +22 dB microwave amplifier, and finally de- 
tected by a spectrum analyzer with an upper frequency 
limit of 46 GHz (Rohde & Schwarz FSU46). The actual 
RF current at the STO is calculated by taking into ac- 
count losses and reflections due to impedance mismatch 
in the transmission line. Losses in our transmission line 
and circulator are characterized by injecting an input sig- 
nal with the microwave source and measuring the output 
with the spectrum analyzer. The reflection at the STO 
is measured with a vector network analyzer and is shown 
in the inset of Fig. [ljb). The scattering matrix element 
Su shown in the figure is proportional to the amount of 
reflection at the STO, which is as high as 70 — 80 % over 
the entire measured frequency range, 0.01-26 GHz. All 
other components in the transmission line, which have 
nominal 50 f2 impedance, give a relative negligible con- 
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FIG. 2: (color online) Frequency and integrated power of 
the free-running STO around the dc bias current values of 
(a,d) 28 mA, (b,e) 31 mA and (c,f) 38 mA. The correspond- 
ing 4 t,l -order polynomial fits to frequency and the 3 rd -order 
polynomial fits to power are shown in solid red lines. 



tribution to the total amount of reflected signal. The 
signal detected at the spectrum analyzer is finally cor- 
rected for the standing waves in the transmission line. 
All data shown in this work have been corrected in order 
to compensate for all these effects. 

The measurements are performed in a magnetic field 
of 10 kOe applied at an angle of 70° to the film plane to 
ensure that (i) the STO operates around its maximum 
output power—, and (m) only the so-called propagating 
modei£~— is excited. This mode has a higher frequency 
than the ferromagnetic resonance mode and shows a blue- 
shift with bias current as confirmed in Fig. [TJa). Fig- 
ure [1] also shows that both the operating frequency and 
the integrated output power (which is proportional to 
the actual precession amplitude of the STO) [Fig. [T[b)] 
are strongly non-linear functions of the dc bias current. 
This behavior is likely related to the excitation of closely 
spaced discrete dynamic modes as the bias current is in- 
creased*^— 

To test different non-linear modulation theories, we 
have chosen to focus on three principally different non- 
linear situations described by three different values of 
d 2 f/dl^ c : zero, positive, and negative, corresponding to 
a drive current of 28, 31, and 38 mA, respectively. These 
three operating points arc shown as dotted lines in Fig. [T] 
The non-linearity can be more clearly seen in Fig. [2] 
which shows the frequency and integrated power of the 
free-running STO around these dc bias current values in 
a range equal to the maximum modulation current. The 
shape of frequency vs current at 28 mA is almost lin- 
ear while it is convex for 31 mA and concave for 38 mA. 
The amplitude sensitivity is also clearly different at these 
current values, as seen from the corresponding plots of 
integrated power in Figs. |3fd)-2(f). Around these oper- 
ating points we modulate the STO using a 100 MHz RF 
signal swept from to 3 mA. The corresponding spectra 
are shown in Fig. [3] as a function of the modulating cur- 
rent amplitude. In all three cases, the number of side- 
bands increases with increasing modulation amplitude. 
In the case of a linear frequency dependence (28 mA, 
d 2 f/dl% c = 0) the carrier and sideband frequencies are 
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FIG. 3: (color online) Frequency modulation (f m = 100 MHz) 
of the STO showing the progressive development of sidebands 
with increasing modulating amplitude I m at dc bias current 
values of 28 mA, 31 mA and 38 mA. Power is expressed in dB 
over the noise floor. The white lines show the calculated fre- 
quency of the carrier and the first-order sidebands according 
to the combined NFAM theory. 



entirely independent of the modulating current (up to a 
modulation current of 2 mA). In contrast, both the car- 
rier and the sideband frequencies show a clear blue-shift 
at 31 mA and a clear red-shift at 38 mA as expected from 
the finite d 2 f/dl 2 c with opposite signs. 

In Fig. SI we show the detailed modulation current 
dependence of the carrier and the first-order sideband 
power with calculated results as described in the follow- 
ing paragraph. While the evolution of the carrier power 
with modulation current does not seem to be affected by 
the non-linearity, both the upper and lower sidebands are 
strongly affected by the sign and the value of d 2 f/dlj c : 
the lower sideband gets markedly stronger than the upper 
sideband for d 2 f/dl% c > (31 mA), and weaker than the 
upper sideband for d 2 f/dl 2 c < (38 mA). The position 
of the maximum sideband power is also shifted up / down 
for the upper/lower sideband. It is noteworthy that this 
shift only depends on the magnitude of d 2 f/dl% c and does 
not change sign when d 2 f/dl 2 c goes from positive to neg- 
ative. Even for the linear case (28 mA, d 2 f/dl% c = 0), 
the power of the two sidebands are unequal. The up- 
per sideband has higher power than the lower sideband, 
as expected from the positive slope of amplitude versus 
bias current in Fig. [2jd) . This case of linear frequency 
modulation provides a strong experimental evidence that 
amplitude modulation is also taking place. 

In order to interpret the observed behavior and esti- 
mate the importance of both the frequency and ampli- 
tude non-linearities, we consider three qualitatively dif- 
ferent models describing (i) LFM, (ii) NFM, and (in) 
NFAM. The latter model is adapted frorn^ and specif- 
ically takes into account non-linearities in both output 
frequency and amplitude as a function of the input bias 
current. 

Since LFM and NFM models have already been de- 
scribed in Re f 14 i 15 , we focus on the details of the NFAM 
model used in our analysis. The instantaneous frequency 




FIG. 4: (color online) Integrated power of the carrier (black 
triangles), and the first-order upper (blue squares) and lower 
(red circles) sidebands for the three different dc bias currents: 
row (a)-(c) 28 mA, row (d)-(f) 31 mA and row (g)-(h) 38 mA. 
First, second and third column shows the corresponding cal- 
culated integrated power (solid lines) as predicted by LFM, 
NFM and NFAM, respectively. The mean square error, % 2 
between the experiment and calculated results of the two side- 
bands improved significantly for NFAM. 



is assumed to depend nonlinearly on the modulating sig- 
nal: 



fi(t) = fc + fcim(i) + k 2 m(t) 2 + k 3 m(tf + 



(1) 



where, m(t), is the modulating signal and the coefficients 
ki represent the z-th order frequency sensitivity coeffi- 
cients. Similarly, the output amplitude, A c is given by 



A c (t) = A + Aim(i) + \ 2 m{t) 2 + X 3 m(t) 3 + 



(2) 



where A^ is ith order amplitude sensitivity coefficient. 
The coefficients ki and Xi are given by the non-linear cur- 
rent dependence of / and A of the free running STO. We 
use sine wave modulation, m(t) = I m sin(2irf m t)i where 
/,„ is the amplitude and f m is the frequency of modu- 
lating signal. The resulting NFAM spectrum becomes^ 



[<*(/ - fc - (« + 2 ™ + 3p + 4g + h)f m ) 
+ 6(f - f c - (n + 2m + 3p + Aq - h)f m ) 
+ 6{f + f c - (n + 2m + 3p + Aq + h)f m ) 
+ S(J + f c - (n + 2m + 3p + Aq - h)f m )} 

(3) 

where ft = k r I m / f m + 3k 3 I? n /4f m , = k 2 I% l /4f m + 
/c 4 /*/4/ m , /3 3 = k 3 lf n /12f m , and f3 4 = fc 4 /^/32/ m are 
frequency modulation indices of different order. 70 = 
Ao + A 2 /^/2, 7l = Ai/ m + 3A 3 / I 3 n /4, 72 = A 2 J^/2, and 
73 = A3-Z^/4 arc amplitude modulation parameters. In 
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TABLE I: Modulation sensitivity coefficients found from polynomial fits of frequency and amplitude of the free running STO. 



Current 


fco 


fei 


k 2 
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Ao 


Ai 
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(mA) 


(GHz) 
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(MHz/ 


(MHz/ 


(pW 1 / 2 ) 


(pW 1 ' 2 / 


(pW 1 / 2 / 


(pW 1 / 2 / 






mA) 


mA 2 ) 


mA 3 ) 


mA 4 ) 




mA) 


mA 2 ) 


mA 3 ) 


28 


20.185 


117 ±1 


1±1 


2 ±0.2 


8± 1 


10.4 ±0.5 


0.9 ±0.07 


-0.2 ±0.02 


-0.03 ± 0.01 


31 


20.545 


147 ±1 


20 ±1 


0.8 ± 0.1 


-1±0.1 


10.9 ±0.6 


-0.5 ±0.07 


-0.15 ±0.02 


0.02 ±0.01 


38 


21.779 


115 ±1 


-22.5 ±0.6 


-3.3 ±0.1 


1.6 ±0.1 


10.8 ± 1 


1.3 ±0.07 


-0.1 ±0.02 


-0.12 ± 0.01 



the above we assumed that the frequency in Eq. (JTJ) is 
non-linear up to fourth order and the amplitude in Eq. ([2]) 
is non-linear up to third order, which is found sufficient to 
describe the experimental data. The frequency spectrum 
S(f) consists of a shifted carrier frequency 

Si = k + k 2 ll + 3fc 4 /4/8 + ... (4) 

and an infinite number of sidebands symmetrically lo- 
cated at ± lf m , where I = n + 2m + 3p + 4q ± h is 
a positive integer identifying the sideband order. The 
NFAM carrier shift is identical to that obtained from an 
NFM model since effects due to amplitude modulation do 
not enter in Eq. (j4|). This shift can be readily calculated 
by means of the polynomial fitting procedure shown in 
Fig. [2j The comparison with the experimentally obtained 
values reveals a good agreement, as shown in Fig. [3] The 
sideband power, on the other hand, is strongly affected 
by the amplitude modulation, through the coefficients 
7t, and can be used to compare the NFM and NFAM 
models. In a 6 mA interval around each operating point, 
we expand the frequency dependence into a fourth-order 
Taylor series, and the amplitude dependence into a third- 
order Taylor series as shown in Fig. [2] The coefficients 
along with their standard errors are summarized in Ta- 
ble |TJ Using these coefficients we calculate the sideband 
power expected from NFM and NFAM, respectively (sec- 
ond and third column in Fig. |4|) and also compare with 
LFM theory (first column in Fig.HJ). 

LFM theory completely fails to describe the strong 
asymmetry between the upper and lower sidebands in all 
cases. In the linear case of 28 mA [Figs. Eta)- BJc)] both 
NFM and LFM theory predict nearly the same behavior 
with equal sideband power since only k\ is significant and 
&2 ~ 0. In contrast, the NFAM model correctly produces 
both the upper and lower sideband power, implying a 
much better agreement, mostly captured by the ampli- 
tude modulation sensitivity coefficient Ai. In fact, the 
mean square error, \ 2 between the experiment and cal- 
culated results of the two sidebands decreases by about 
80 % for NFAM theory compared to LFM. In the two 
non-linear cases, the NFM model captures the change 
in sign of the sideband asymmetry, given by the sign 
change in hi , but only yields a partial improvement com- 
pared to LFM. On the contrary, when the amplitude sen- 
sitivity coefficients are also taken into account the agree- 
ment of the calculations with experiment is essentially 



perfect. This agreement is only obtained when both fre- 
quency and amplitude non-linearities are accounted for; 
both &2 and Ai are significant. For 31 mA (38 mA), 
the mean square error between the experiment and cal- 
culated results of the two sidebands decreases by about 
85 % (83 %) for NFAM theory compared to LFM and 
about 10 % (36 %) compared to NFM theory. We empha- 
size that none of the presented calculations involve any 
free parameters and are completely determined by the 
experimentally measured nonlinear current dependences 
of the free- running STO. The agreement with NFAM was 
also found to be valid for a range of lower modulation fre- 
quencies (down to 40 MHz) over the entire range of dc 
bias currents. Thus our results show that, as long as both 
non-linearities are accounted for, the proposed scheme of 
combined modulation is able to accurately predict the re- 
sulting sideband powers and frequency shifts over a wide 
range of varying operating conditions. Consequently, the 
STO behaves as an ordinary RF oscillator and should 
lend itself to communication applications. 

In conclusion, we have carried out a detailed modu- 
lation study on a nano-contact STO. In particular, we 
have studied the impact of different levels of frequency 
non-linearity. In the non-linear cases, both carrier and 
sidebands frequencies are shifted as a function of the 
modulation current. Both frequency and amplitude non- 
linearities produce a significant asymmetry in the power 
of the upper and lower sidebands. We find that a com- 
bined non-linear frequency and amplitude modulation 
model can accurately describe all our experimental data 
without any adjustable parameters. The modulation of 
an STO is therefore predictable and independent of the 
complex mechanism behind the non-linearity. The results 
are significant for the continued development of commu- 
nication and signal processing applications of spin torque 
oscillators. 
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